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Slit shaped apertures are usually used in small-angle x-ray scattering measurements in order to 
obtain easily measured intensities of scattered radiation. As a result, the scattering intensity at a 
given angle determined by the camera, the center of the sample, and the central incident x-ray beam 
is not simply related to the scattering from the sample only at that angle. The experimentally deter- 
mined intensities I{x) are related to the true scattering intensities by the following integral equation: 



Hx)-- 



i[{x' + tyi^-W{t)dt. 



This integral equation has been previously solved only for certain simplified functional forms for JV{t) 
and I{x). In this paper, a formal procedure is developed for calculating I{x) from the observed angular 
meas-urements, which does not necessitate making any a priori assumption about the form of JV{t) 
and I{x). 
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1. Introduction 

Narrow slits are often used in small-angle x-ray 
scattering measurements to collimate the incident 
beam sufficiently to allow scattering measurements 
extremely close to the incident beam. Compared to 
a pinhole of diameter approximately equal to the slit 
width, the slit system increases the available scatter- 
ing energy. However, the scattered radiation meas- 
ured in a given solid angle can no longer be simply 
described as the scattering coming from a point source. 
Instead, the actual angular scattering measurements, 
which are summations of the scattering contributed by 
each scattering center along the height of the sample 
slit, must be resolved to give the scattering of an infi- 



nitely small area at the center of the sample slit. In 
this paper, a formal procedure is described for calculat- 
ing the true scattering cross sections from the observed 
angular measurements when long but very narrow- 
width slits are used for collimation. 

A diagram for the slit system of a small angle x-ray 
system is given in figure 1. This is a modification 
of the diagram given by ShuU and Roess [1].^ The 
plane of the quadrangle txyr is perpendicular to the 
camera axis which is in the vertical plane containing 
the collimating slit and the sample slit. The col- 
limating slit nearest the x-ray tube has length 2/], the 
effective length of the slit at the sample is 24, and the 



' Figures in brackets indicate the literature references at the end of this paper. 




Figure 1. Diagram of a slit collimation system. 
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scanning slit at the detector is of length 2yo. The 
scanning slit and sample slit are shown in figure 1 to 
be in a plane making an angle with the plane con- 
taining the camera axis and the coUimating slit. Let 
us consider a ray which makes an angle 7) with the 
camera axis and after it hits the sample is scattered 
by an angle e with its original direction. Since the 
angles e, 6, and r) are in practice small, they can be 
replaced by the following ratios: 



e=xiL, 

and 7] = (t — s)IL, 

Since r^ = x^ -\- {t — yf , it follows that 



e2 = ^2_^ 



7) + 



s — y 



(1) 



(2) 



Let I{6) be the intensity striking the scanning slit cen- 
tered at ^, and let I{e)dr] represent the increment in 
intensity of scattering in an element dr\. The fol- 
lowing equation defines I{S) by averaging /(e) over the 
entire scanning area. 

1 r^o ^2 ris+ii)iL; 
m = W- dy ds I{e)dr). (3) 

Equation (3) is the one obtained by Shull and Roess []],. 
The triple integral in eq (3) can be reduced to a 
single integral. One way to do it is by the method 
employed by Schmidt [2]. By changing y, 5, and r) 
into the dimensionless variables 



yi=ylLu si=s f^+^j, ^71 = ^7- 



s 

u 



the hmits of the integrations are replaced by constants 
which are descriptive of the geometrical parameters 
of the slit system. One can then integrate the re- 
sulting expression over 5i and y\. After substituting 

ti = -j- for 171 +5i, the resulting integral equation is 



m 



Jo 



/[(6»2 + /;f)'/2]r(f,Mfi. 



(4) 



However, W{t\) is a comphcated function of the geo- 
metrical parameters of the system. In general, W{t\) 
is a quadratic function of tu represented by different 
functions for different intervals of ti. Also, W{ti) 
must be premultipHed by a unit step function which 
becomes zero when 



double integrals of the general form 

I I g{'r]\-si-yx)dsidyu 

with the ranges of integration being determined by 
the parameters of the system. In other words, W{t\) 
in eq (4) results from the projections of the variables 
5i and y\ onto the plane of registration, corrected for 
a possible nonuniformity in the distribution of intensity 
of scattered ratiation along 7)\. The only restrictions 
imposed are on the variables 5i, 171, and ji. These 
angles must be sufficiently small to justify their re- 
placements by ratios of coordinates (eq (1)). 

In the process of determining the formal solution of 
the integral eq (4), one has to calculate derivatives of 
W{t), and for this reason, it would be more advan- 
tageous to start with a continuous expression for this 
function. This difficulty is eliminated if one starts 
with the experimentally determined W{t) or by using 
the smoothed curve obtained from the calculated 
segments. 

Starting with eq (3), but omitting the slit coordinate 
171, we have the following equation, with w{t) repre- 
senting the intensity distribution in the t direction 



I{x) 



J — oc J — 



wit)I{[x'^it-yryi^dtdy. (5) 



Interchanging the order of integration and making use 
of the fact that W(t) = W{— t), we obtain 



Jo 



with 



J- 



w{t — y)dy. 



(6) 



(7) 



W{t) can be determined by measuring the intensity 
in the plane of registration. This is done by taking 
various intervals equal to 2yo along the t axis. These 
intervals are equal to the length of the scanning slit. 
The finite length of the scanning slit restricts the range 
of the variable y in eq (7). Therefore, a unit step 
function h{y) has to be introduced, with h(y) = for 
y>k and h{y)=l for 0<y<A:. This process will 
modify Wit) into a convolution integral of w{t) and h{t), 



i 



rU)= uiy)h{t-y)dy. 



(8) 



t> 



yo+h+{li-\-l2)-r 



The advantage of starting with eq (3) is that £3171 is 
the coordinate of the length of the slit nearest to the 
x-ray tube. Thus, in the more general case when the 
intensity along 171 is not uniform, but is given by a 
general distribution function g(y]i), eq (4) will still be 
obtained. However, W{t^ will be represented by 



The only known solutions of the integral eq (6) were 
obtained for certain special forms of W{t) and of 

l{x). Guinier [3], Fournet [4], and DuMond [5] have 
solved this integral equation for shts of infinite height 
{W{t) = 1). Their method was extended by Kratky, 
Porod, and Kahovec [6], and by Ceroid [7] to finite 
beams where W{t) could be represented by a step 
function, or a Gaussian distribution, or a particular 
trapezoidal distribution. Shull and Roess [I] have 
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presented a solution for an arbitrary W{t) if the scat- 
tering curve I(x) can be resolved into a packet of Gaus- 
sian functions. The earlier attempts to solve eq (6) 
have been recently reviewed by Beeman, Kaesberg, 
Anderegg, and Webb [8], and by Kratky, Porod, and 
Skala[9]. 

In the following section a formal solution to eq (6) 
is presented which does not require any a priori as- 
sumptions about the particular forms of I{x) and of 



2. Solution of the Integral Equation 

2.1. General Considerations and Discussions 

Equation (6) is rewritten as 

yI{y)W[if-xyi']dy (9) 



Explicit Solution: 



.H2/; 



I(x) = 



if-xjl^ 



To solve this integral equation in one of the standard 
forms it is essential to remove its singularity, since 
the kernel of the equation is infinite at y = x. The 
transformation of this equation, which removes this 
singularity foHows the method employed in connection 
with the celebrated Abel's integral equation. A new 
variable in a suitable form is introduced, and both 
sides of the resulting integral equation are integrated 
over the variable x between the value of the new vari- 
able and infinity. This, combined with the inter- 
change in the order of integration of the right-hand 
side of the resulting integral equation, through the 
use of the Dirichlet formula, has the effect of remov- 
ing the singularity. The new kernel is finite every- 
where, but it is replaced by a bounded integral whose 
integrand involves the original distribution function 

The solution thus obtained is an implicit one, since 
the left-hand side of the integral equation involves 
the integrated form of the experimental function I(x). 
In order to obtain an explicit form for the integral 
equation, two additional steps are required: (1) Solution 
of a subsidiary integral equation for another function, 
which essentially compensates for the sht weighting 
function JF, (see following section), and (2) differentia- 
tion of the experimental intensity function lix) with 
respect to x. The details of the solution are given in 
the following section. Here we give the two solutions, 
the imphcit and the exphcit one in their final forms: 
Implicit Solution: 



f 



/. 



I[{v^ + uyi^dv = 2 yI{y)K[{f-u^yl^dy, (10) 



with 



K( 



^»=ra 



TT jo 



(y-h/i2)i/2 



F{y)dy. 



(11) 



The function F satisfies the following Volterra equa- 
tion of the second kind: 



F{y)=l--yf " F{z)K\{f - z^l'/^Jfiz, 

TT Jo 



(12) 



K{i^) 



-I' 

I' Jo 



xW'{vx)dx 



2\l/2 



'-=^f-z'. (12a) 



The kernels of the imphcit and expHcit solutions, eqs 
(10a) and (12a), are similar in their general forms, 
except that the first one involves the slit correction 
function W, and the second one involves its derivative. 
The implicit solution is suitable only with the use of 
a digital computer. Both sides of eq (10) are solved 
simultaneously as a set of simultaneous linear equa- 
tions through a matrix inversion method. This method 
can be employed also in order to obtain F{u) from 
eq (12). However, the latter equation can be solved 
by expansion of F{u) in a Liouville-Neumann series. 
These series will converge rapidly, if a good initial 
assumption is made about the approximate form of 
the function F{u). 

Summarizing, one can demonstrate that the solu- 
tion to the problem of slit correction always leads to 
a single integral equation provided that all angular 
variables of the x-ray camera-slit system are small. 
The solution of this integral equation invariably leads 
to a solution of the Volterra equation of the second 
kind which has for its inhomogeneous part a constant 
term, and whose kernel involves an integrated expres- 
sion of the derivative of the experimentally determin- 
able slit correction function. 

At the present time a computer program for the 
numerical solution of the eqs (10-12) is being 
developed. 

2.2. Derivation of Equations (10-12) 

In Abel's integral equation, the kernel K{x, y) is of 

the form — , with < a < 1 and x> r. For 

{x — yY ^ 

a— 1/2, the singularity of the kernel is removed by 
multiplication of both sides of the integral equation by 

; TTTT and integrating with respect to x. This is 

{u — xY''^ 

done to make use of the fact that 



W{vx) 

V-2V/2 



1^ dx 

J, [(x-y){a-x)VI^''' 



(13) 



dx', 1^2 y2 ^2 ^Q^^ ^^^^ r[(/-Jc2)i/2]-/C(jc, y), eq (6) is written in 
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230-736 0-66— 2 



Abelian form: 



/( 



^>=^r 



yI(y)K{x, y) 



dy. 



(14) 



For this particular equation, the transformation which 
removes the singularity is based not on eq (13), but 
on an equivalent form: Use is made of the fact, that 



I: 



X dx 



u [{y^-x^){x^-u')^ 



1/2 



77 

"r 



since eq (15) can be shown to reduce to eq (13). 
therefore multiply both sides of eq (14) by 

xdx 



(15) 
We 



2 _ „2\l/2 



(x 



and integrate over x from u to <^, obtaining 



r^ xi{x)dx, r°° r^ 
ju (x'-u'yi'~^ ju J. 



xyI{y)W[{y^-x^yl^dydx 
[{y^-x^){x^-u^W 



Ju 



yI{y)Uu, y)dy, (16) 



nth 



, fv xW\y^-x^yi^]dx 
^'^"'''^"Ju [{x^-u^){y'-x^Wl^' ^^^^ 

In obtaining eq (16) and (17), the order of integration 
was exchanged using the following Dirichlet formula: 

Too Too Too fy 

dx G{y, x)dy= dy G{y, x)dx. (18) 

JU Jx Ju Ju 

The implicit solution (eq (10-lOa)) is obtained from 
(16) and (17) in a straightforward way: we replace 
y^ — x^ by t^. The kernel (17) simplifies to 



u 



u,y) = l 



\y^-un' 



W{t)dt 



(3^ _ ^2 _ ^2)1/2' 



(19) 



and eq (10-lOa) follow directly. To obtain the ex- 
plicit solution of eq (14), we multiply the integrands of 
both sides of eq (16) by F[{x^ — u^y>^\ F{z) is an un- 
known function to be determined subsequently. In 
the right-hand side of eq (16) K2{u^ y) replaces Ki{u^ v) 
and is given by eq (20): 



{y xFlix^-u'^yi^Wiiy'-x^y 

'^ '^' L [iy^-x^){x^-u')y^ 



l^]dx 



(20) 



Both sides of the resuhing eq (16) are now differ- 
entiated with respect to u: 



d r°° xJ(x 
duju ( 



dx 



)F[{x^-u^yi^] 
(x2-u2)i/2 

= -2uI{u)K-2{u, u) + 2 



f\/(r)^^^.y. 

J u 



du 



(21) 



We need a function F(z), {z^ = x^ — u^), which leads to 
K2{u, y) = X, independent of u and y. If .^i(x, y) is 
constant, F{z) then would be a constant also, equal 
to 2/77. We would then have an Abel-related integral 
quation with a well-known solution for I{x). With 
K'ziu, y) = X, the solution of eq (21) is readily obtained: 



Ku) 



2Xjo 



i[iz^-+u^y'^] 

2\Jo (z2-hf/2)l/2 



F{z)dz. 



(23) 



Equation (23) was obtained by Kratky, Porod, and 
Kahovec [6], using a different method of derivation. 
At first glance, the difficulty we encountered in 
solving the original integral eq (14), which is brought 
on by the singularity of its kernel, does not seem to be 
removed. The integral equation for the new function 
F{z) still possesses a similar singularity. However, 
the integral equation for F{z) (eq (20)) is considerably 
simpler than the original integral equation for I{x) if 
K'ziu, y) = constant. The same procedures which 
were used to remove the singularity in the kernel of 
eq (14) lead to a straight-forward solution for F{z) by 
reducing it to a Volterra-type equation of the second 
kind. Equation (20) is rewritten as 



ny2-M2)i/: 

f 



F{z)W[{f-u^-z^yi^] 



if- 



,2 _ ^2M/2 



dz = \. 



F{z)W{{s''-2:'yi^] 

(52-z2)l/2 



dz = X, 



(24) 



with s^ = y^ — a^. 

Both sides of eq (24) are then multipHed by 
s{i^ — s^Y^^^ and integrated over s from o to u. The 
following equation is obtained: 



17' 

Jo jo 



sF{z)W[{s''-zJl^dzds 

[(5--Z^)(«2-,2)],/2 



Jo W 



sds 



.^2)1/2 



-Xu. 
(25) 



We again interchange the order of integration in the 
l.h.s of eq (25), using the Dirichlet formula to obtain 



|V(z)rf2jj 



sWlJs^-z^yl^ds 
[{s^-z'){u^-s^)] 



1/2 



ku. 



(26) 



Next we subtract and add lFiz = s) = W{0) and inte- 
grate eq (26) to obtain 



■r(0) 



( F{z)dz - r F{z)M{u, z)dz = \u, (27) 
Jo Jo 



Mi 



u,z)= I 



s{WiO)-W[{s^-z^y'^]} 

[(52-z2)(a2-52)]l/2 



ds. (27a) 
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After differentiating (27) with respect to u we obtain 
the following VoUerra-type equation of the second kind 

F{u)= Tw.^. + jw.^x — 1 F(z)rfz, (28) 

since M{u, u) = 0. 

Equation (28) can also be further simpHfied. If we 
set r(0)-l and \ = 7t/2, then 

F(w) -= 1 + - 7^—^ F{z)dz. 29) 

77 Jo du 

Equation (27a) is integrated by parts, 



ith 



Kiu.z) = —^^ 

(u^-z-) Jo 



^2)1/2 



[iu 



yff^'{y)dy 



(32) 



The kernel of eq (31) is finite for all values of z, includ- 
ing z= u. 

Equations (12-12a) follow directly by exchanging 
variables, ii being replaced by y. 



3. References 



'TT , 



M(il, z) = ^[l-r[(^2_^2)l/2] 



+ 



'^ sW'is'-z^yi''] 



is' 



^2X1/2 



Therefore, 
dM{ 



s2_^2yl2 



ds. 






du 

Integral eq (29) assumes the following form 
2 



F{u)=} 



rr Jo 



F{z)K{u. z)dz. 



(31) 
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